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ABSTRACT 

In this paper we present and discuss a new Monte Carlo approach aimed at correcting 
the observed Supernova (SN) rates for the effects of host galaxy dust extinction. The 
problem is addressed in a general way and the model includes SN position distribu- 
tions, SN light curve and spectral library, dust properties and distribution as input 
ingredients. Even though the recipe we propose is in principle applicable to all SN 
types, in this paper we illustrate the use of our model only for Type la. These rep- 
resent in fact the simplest test case, basically due to their spectroscopic homogeneity 
which, to a first approximation, allows one to treat them all in the same way. This 
test case shows that the final results do not depend critically on the spiral arm dust 
geometry, while the total amount of dust, its properties and the size of the galactic 
bulge do have a strong effect. With the availability of more complete spectral libraries 
and a more accurate knowledge of SN spatial distribution, the method we propose 
here can be easily extended to Core Collapses events. 
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1 INTRODUCTION 

Until a few years ago, observational measurement of super- 
nova (SN) rates were scanty. The only available measure- 
ment s were those for local SN rates dCappellaro et allll997l 
I1999T) whilst nothing was known about the evolution with 
look-back time. The turning point was represented in the 
late nineties by the demonstration that Type la SNe can be 
used as precise distance indicators, p articularly well suited 
to probe the ge o metry of Universe dSchmidt et alJ Il998t 
iRiess et al.lll998l : IPerlmutter et al.lll999l) . The striking re- 
sult that the Universe is accelerating has indeed stimulated 
the kick off of several new SN search projects, aimed at us- 
ing Type la SNe as cosmological probes both to improve the 
determinations of the density parameters (Qm, ^a) as well 
as to trace the evolut ion of the expansi on factor through the 
Universe history fsee lLeibundeu J200ll for a review) . Before 
that time, studies of SN rates and their evolution with red- 
shift were mainly theoretical speculations. Today, the few 
published measurement are based o n data from SN searches 
desig ned for cosmological purposes jPain et al.lll99rl 120021: 
iTonrv et al.l2003l:lDahlen et al.l2004l). The only exception is 
represented bv lCappellaro et alJ (|2004), a unique SN search 
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specifically designed to measure the SN rate at intermediate 
redshifts and with no bias towards la's. 

A major point of concern when deriving observational 
estimates of SN rates is the correction for all possible bi- 
ases affecting the observed number counts. Nowadays, the 
main source of bias is caused by the effect of dust extinc- 
tion in the SN host galaxies. This is indeed severe, partic- 
ularly for core collapse (CC) SNe, which are preferentially 
located in dust-rich environments. On the contrary, the ef- 
fect for Type la SNe is usually assumed to be milder, par- 
ticularly for those exploding in elliptical and early spiral 
galaxies, which are thought to be almost dust-free environ- 
ments (see e.g. IRiess et alJll996T) . Nevertheless, there are 
several examples of Type la SNe spanning a large range 
of extinction fro m moderate to high values: e.g. SN 1994D 
A v = 0.186 mag JPatat et al.lll99rj) . SN 2002er A v = 0.62 
mag dPignata et alJl2004h . SN 2002bo A v = 1.364 mag 
jBenetti et al J 1200411 . SN 2003cg Ay ~ 2.3 mag (Elias De 
La R osa, private communi cation) and SN 2002cv Av — 8.0 
mag <Di Paola et al.ll2002l) . There are several methods to 
estimate the amount of extinction suffered by a SN: they 
require either the availability of the SN light curve (LC) 
in at least two passbands (around maximum light or be- 
tween 30 and 90 days past B maxim um light) or at least 
one spectrum with good S/N (e.g. see lPignata et al.ll2004l . 
for an application). These techniques are nevertheless use- 
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less when one needs to estimate the bias introduced by host 
galaxy extinction on a sample of SNe discovered by a search 
program. Indeed in such a case, the valuable missing piece 
of information is a statistical correction factor that would 
include, in the SN number counts, the contribution of those 
objects that were extinguished by the dust of their parent 
galaxy and consequently weakened below the search detec- 
tion limit. 

So far the only work aimed at studying the effect 
of extinction on the radial distributions of SN proper- 
ties in their parent galaxy is rep resented by the model of 
lHatano. Branch fc DeatorJ ljl99Sl. hereafter HBD98). Using 
a Monte Carlo approach based on assumed distributions of 
dust and SN populations, they derived the expected extinc- 
tion distributions and projected radial distributions of both 
Type la and CC SNe. In particular, they adopted a sim- 
ple description of the dust distribution from which they de- 
rived the total £?-band extinction per kpc in the plane of 
the disk j a(r), which is a linear function of r (see Eq. 1 
and 2 of lHatano et al.lll998t) . In such a parametrisation, the 
dust density rises from zero in the centre of the model to 
a maximum value at 5 kpc and then falls to zero at 17.5 
kpc. They found that as the galaxy inclination increases, 
the projected distribution's peak moves toward the central 
regions. Moreover, they conclude that the effects of extinc- 
tion induced by the host galaxy are likely to be the cause 
of the observed difference in absolute magnitude dispersion 
of Type la SNe within and b eyond 7.5 kpc, confirming the 
findings of IWang et all Jl997t) . 

Despite the obvious interest in the observed spatial de- 
pendence of SNe properties, in this work we will rather focus 
on deriving a recipe to correct the observed SN rate for host 
galaxy extinction. The HBD98 model has been used to de- 
rive e xtinction corr ections for both loc al (ICappellaro et all 
119991) and high-z flbnrv et all 120031) SN rates. In par- 
ticular, ICappellaro et all il999t) have carefully tested the 
HBD98 approach to derive the extinction corrections as a 
function of galaxy inclination. Interestingly enough, they 
found that the model appears to overcorrect the SN rate 
of distant galaxies suggesting that a more detailed mod- 
elling of the dust distribution and properties is possibly 
required. All modern measurement of SN rates at high-z 
currently a vailable have either not been corrected for dust 
extinction (ICappellaro et all Eool iPain et alJ Il996t l2002t 
iTonrv et al.l2003l) or were corrected jDahlen et al.l2004l) us- 
ing the HBD98 approach. 

For these reasons, and given the importance of this 
topic in modern astrophysics, we have decided to address 
the problem with a more detailed approach. Our final goal 
is to construct a reliable model which is able to provide 
the "bias factor" needed to correct the observed SN rates 
for dust extinction. In this framework, we started from the 
Monte Carlo approach of HBD98, which clearly has proven 
to be effective, and tried to construct a more general model 
able to cope with the variety of dust properties and distri- 
butions that are observed in different spiral galaxies. This is 
indeed a key feature to enable a reliable determination of the 
dust bias factor for a real SN search, where the galaxy sam- 
ple is clea r ly a m ixture of types and morphologies. Recently 
IComminsl (l2004h developed some models for extinction of 
Type la SNe due to dust in spiral galaxies. The main pur- 
pose of his work is attempting to account for observational 



effects that play a significant role in searches for high red- 
shift Type la SNe. In particular, that paper investigates the 
expected dependence of the observed amount of extinction 
in high redshift Type la SNe rather than trying to quantify 
the effect of dust extinction on the measured SN rate s. In 
this respect, we feel that the study of IComminsl <l2004l) and 
the present work are complementary. Nevertheless we will 
discuss the main differences between our model and Com- 
mins's one in Section |21 

In this paper we present the model and the general 
recipe to derive the extinction correction factor while the 
application to the computation of extinction-corrected SN 
rates will be discussed in a companion paper. In Section 
we outline the model and the assumed spatial distributions 
of dust and SNe. Section |H] presents model predictions of 
the extinction distributions and projected radial distribu- 
tions together with their dependence on model parameters. 
As an illustration we show the effects of dust for a single 
test case in Section 2] Section |S] presents a discussion of the 
physical parameters which have a major impact on the ex- 
tinction correction factor. Finally we summarise our results 
in Section [S] 



2 METHOD 

The method we follow in this work is based on the Monte 
Carlo technique and it consists of a few basic steps grouped 
in two main stages: the Monte Carlo and the Integration 
part. In the first stage Type la SN positions inside host 
galaxies are generated according to different distribution 
functions (see Sec. 12. H . After that, random SN phases and 
lines of sight are uniformly and isotropically generated, re- 
spectively. These steps complete the Monte Carlo stage of 
the simulation. The Integration stage uses the Monte Carlo 
data to compute the SN magnitudes. First the total col- 
umn density along the line of sight associated to each object 
is computed by integrating the dust density function (see 
Sec. 12.21 along the specific line-of-sight. Finally, the appar- 
ent magnitudes of SNe are computed integrating template 
SN spectra, redshifted to the suitable redshift z, over the 
selected passband V in the presence, or a bsence, of dus t ab- 
sorption using Draine's extinction laws JPrainell2003l (see 
Sec. 12.31 . In this way we can study how the distribution of 
SN extinction depends on different geometries (e.g. galaxy 
inclination, different SN distributions inside the host galaxy, 
etc.) and different physical parameters (e.g. face-on cen- 
tral optical depth, etc.). In this work we have always used 
Draine's extinction law for Rv = 3.1 but the model imple- 
mentation considers also the cases Rv = 4.0 and Rv = 5.5. 
In principle, any another extinction law for which the extinc- 
tion coefficients as a function of wavelength are provided is 
suitable as well. We note that, from the analysis of individ- 
ual SNe, in some cases values of Rv have been found as low 
as 2.0 (Elias de la Rosa, private communication). Though 
the approach we used is general, we decided to limit our 
discussion only to Type la SNe both because of the lack of 
enough data (LC and spectra) for each CC SNe type (lb, Ic, 
IIP, IIL, etc.) and because of the present great interest in 
Type la SNe due to their cosmological applications. 

Throughout this paper we will assume the standard 
cosmological model (Ho = 71 km s _1 Mpc -1 , Q,m ~ 0.3, 
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Qa = 0.7). The adopted values for the other parameters 
will be discussed in the following sections. 

2.1 Type la SNe distribution 

It is generally accepted that Type la SNe originate from the 
thermonuclear explosion of accreting white dwarfs (WD) in 
binary systems. The debate on the nature of the progeni- 
tor systems is still active. In particular there are two main 
scenarios: (1) the single degenerate, in which the progen- 
itor WD accretes hydrogen-rich matter via mass transfer 
from a companion main sequence or low m ass RGB star 
JWhelan fc Ibenlll984l: llben fc Tutukovlll984h ; (2) the dou- 
ble degenerate, in which the binary system is composed of 
two WDs, the exploding WD reaches the Chandrasekhar 
limit and explosive c arbon ignition occurs at its centre 
jThen fc TntukovllT98l. 

It thus seems reasonable to assume that the old galaxy 
stellar component is a good tracer of the Type la SNe spatial 
distribution. The stellar component in a spiral galaxy can be 
described by two structures: a spheroidal bulge and a disk. 
To model these two components, we mainly followed the 
method employed in radiative transfer codes used to derive 
extincti on and polarization effects in dusty spiral galaxies, 
(see e .g. lBianchi. Ferrara fc Giovan ardi 19961. IXilouris et alJ 
1999). We hence used the luminosity density distributions 
for the two galaxy components to describe the spatial num- 
ber density of Type la SNe, which we will denote as pi a (see 
also IComminslEool) . 

To reproduce the variety of spiral Hubble types, we have 
considered different bulge-to-total luminosity ratios (B/T) 
in the range (resembling an Sd galaxy) to 0.5 (resembling 
an Sa). Given the parallelism between the luminosity distri- 
bution and the Type la SN spatial distribution previously in- 
troduced, once the total number, Af* a of simulated Type la 
SNe has been fixed, a fraction B/T is generated according to 
the bulge distribution and the remaining fraction (1 — B/T) 
is generated according to the disk distribution. 

We note that the properties of the assumed spatial dis- 
tribution of Type la SNe are among the most important 
parameters of our model. Although the adopted spatial dis- 
tribution is rather simple, it takes into account the diffent 
scale height of dust and SNe suggested b y the observations 
JPella Valle fc Panagialfl992l: |Patatll2005l) . A more detailed 
treatment of the differences between the distribution of the 
stellar component and Type la SNe could improve the model 
but requires a careful analysis of the differences in the distri- 
butions between single stars and primordial binaries which 
is beyond the scope of the present study. 



2.1.1 



Bulge distribution 

For the bulge, fol lowing [ Bianchi et alJ 
Jaffe distribution (|Jaffelll98!l )~~ 



PlaW = 



pi 



(r/rt) [l+(rM)]- 



, we adopt a 



(1) 



where r is the physical distance from the galactic centre 
and ri a is the scale radius. The volume integral of the Jaffe 
distribution converges for r — > oo but clearly, for numerical 
reasons, we need to introduce a truncation radius r! 



nr Ia . Without this truncation, at large radii the number of 
bulge SNe would indeed dominate once again that of the 
disk. In this work we used rj a = 1.87 for the bulge scale- 
length and n — 8 for the truncation radius. This value of 
r Ia , which corresponds to an equivalent radius R e = 1.60, 



is clo se to the one inferred for our Galaxy jBianchi et alJ 
1996). 



2.1.2 Exponential disk distribution 

As for the bulge, we assume that the spatial distribution 
of stars in the disk is traced by the luminosity density. As 
a fir st approximatio n we can consider a double exponential 
law <Freemadll97ch . namely 



pL(R, 



Po exp 



R 



(2) 



where R and z are the cylindrical coordinates, Rf a is the disk 
scale- length and Zd the disk scale-height. Following the same 
approach as in the bulge case, we introduce a radial, = 
nRf a , and a vertical, = mzd, external truncation. In 

this work we used rf a = 4.0 kpc and zf a = 0.35 kpc which 
are similar to the values inferred for the old disk population 
of the Gala xy. For the outer e dge truncation we adopted 
n = m = 6 jBianchi et al.lll996h . 



2.1.3 Perturbed exponential disk distribution 

As the aim of this work is to improve the calculation of 
the statistical properties of Type la SNe extinction by host 
galaxy dust, it is crucial to simulate their spatial distribution 
inside the gala xy as realistically as possible. In particular, 
in their study, McMilla n fc Ciardullol (I1996T) considered a 
sample of 35 spectroscopically identified Type la SNe and 
measured their positions relative to the spiral arms from 
a homogeneous data set of V and I images. Interestingly 
enough, they found that Type la SNe occur closer to the 
arms than a perfectly randomised disk population. However, 
they note that Type la SNe are not statistically distinguish- 
able from the distribution of stellar light around the arms 
(which should be considered as a tracer of the stellar den- 
sity). 

We decided thus to implement in our model both cases: 
(1) no concentration closer to arms, which we achieve by dis- 
tributing the SNe in a homogeneous exponential disk (Eq. 
[3J; (2) concentration closer to arms, which we achieve by 
simulating a spiral structure as a perturbation of a homo- 
geneous exponential disk. A realistic s piral structure can be 
repre sented by logarithmic spiral arms feinnev fc Merrifieldl 
1998) as a perturbation o n a plain exponential disk. Follow- 
ing Misiriotis et alJ (l2000Tl , for the perturbed distribution we 
adopt the following expression: 



pL(R, 



1 + w sin [ ^ a log R — N a d> 
\ tanp 



(3) 



' m a x 



where R, z and <j> are the cylindrical coordinates. pj a (_R, z) is 
made up of two components: the plain exponential disk pj a , 
given by Eq. |21 and the spiral perturbation represented by 
the term in square brackets. The three parameters govern- 
ing the spiral perturbation are: w, the spiral perturbation 
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amplitude parameter (w — correspond to the plain expo- 
nential disk); p, the pitch angle (low p values produce tightly 
wound arms) and N a , the number of spiral arms. 

Though we have implemented in our model the per- 
turbed spiral distribution described by Eq. |SJ all the re- 
sults discussed in this work are based on simulations in 
which the disk component was described by Eq. |2] We 
believe indeed that, t hough the evidences presented by 
iMcMillan fc Ciardullol |l996) appear quite convincing, ob- 
servational data is still too scanty to permit modelling with- 
out "blind" assumptions. We thus prefer to limit the study 
to the simpler case to reduce the number of low-impact 
model's free parameters. 



2.2 Dust distribution 

The exponential disk is commonly used in radiative trans- 
fer Monte Carlo codes of galaxies because it is easily im- 
plemented numerically and, moreover, it makes simulations 
very efficient in terms of computation time. In our approach 
we do not solve for the radiation transfer and thus we can 
implement more realistic dust distributions including a spi- 
ral arm structure, still keeping an affordable computation 
time. The dust density distribution p~d(R, z, </>) is described 
by an expression analogous to Eq.[3] where the choice of the 
spiral perturbation parameters (N a , p and w) and the plain 
disk scale length and height (Rp and z%, analogous to Rf a 
and zfg, in Eq. |5J are obviously independent from that used 
for Type la SNe distribution. In analogy with the Type la 
SNe case, we have introduced a radial, R^ ax = tiRd, and 
a vertical, Zmax = mzn, external truncation. In this work 
we adopt Rd = 4.0 kpc and zd = 0.14 kpc and the same 
truncation factors as for the Type la SN distribution (see 
Sec. 12.1.21 . The total content of dust in a galaxy is then 
normalised in terms of the optical depth along the galaxy 
symmetry axis: 

/ + oo 
p d D (R = 0, z)dz = p d [2z d e- m C ext {\f )] (4) 
- oo 

where m is the disk vertical external truncation introduced 
before and the subscript V indicates that the normalisation 
depends on the adopted passband, which is characterised by 
the value of the extinction coefficient at its effective wave- 
length A^?. Once r-p(O) has been fixed, Eq. 2] provides the 
normalisation pfj f° r the dust distribution. 

Eq. 0]is rigourously valid for the plain disk only, but it 
has been used also for the perturbed exponential disk be- 
cause, in such case, the normalisation has not an analytical 
expression. 

We note here that HBD98 did not allow the dust den- 
sity to increase to the very centre of the galaxy. Conversely 
our models do not adopt any internal cut-off in the dust dis- 
tribution. It appears that the central ho le is a well know fea- 
ture o f the Milky Way dust distribution jDrimmel fc Spergell 
l200ll) . We note however that usually such a central cut-off is 
not adopted by modern radiative transfer models used to fit 
the optical and near infrared surface brightness of observed 
galaxies. Therefore, for the lack of compelling evidence to 
the contrary and for the sake of generality, we assume a 
smooth dust distribution without a central cut-off as de- 
scribed above. This difference may partially justify some of 



the diversities we found between our results and HBD98's 
ones (see Section f3.H , 

2.3 Type la SN apparent magnitudes 

To compute the SN apparent magnitudes one can follow two 
different approaches, one based on SN light curve and the 
other based on SN spectra. In the first approach one has to 
select a particular SN light curve template, thus choosing a 
particular Amis parameter (defined as the B band magni- 
tude difference between maximum light and 15 days after). 
Other parameters are the absolute magnitude at maximum, 
its intrinsic dispersion and the isT-correction which depends 
on redshift, filter and SN type and phase. Finally one has to 
derive the total extinction along the line of sight and in the 
selected passband. 

In the second approach, which we choose to follow, one 
has to select a set of template spectra, covering a reasonable 
phase range. Then, for a given SN phase, the spectrum is 
extinguished according to the total optical depth along the 
line of sight, and integrated over the passband V to derive 
the apparent flux Fj>: 
p+00 

/ F(X) cxp [-N d Cg t (X)]T v (\) dX 

Fv = ^ (5) 

/ T v (X)dX 
Jo 

where F(X) is the (input) observed SN spectrum, C ( ! xt (A) 
is the extinction cross section per H nucleon (cm 2 /H) as a 
function of wavelength for a given value of Rv (see Sec. l2.5H . 
Tp is the total transmission for the passband V and Nd is the 
column density obtained by computing the integral of the 
dust density function along the line of sight. The apparent 
magnitude of the SN is finally computed converting the flux 
to magnitudes and applying the photometric zero point zpv, 
namely: 

m v = -2.5 {log [F V {1 + z)] - log (zp P )} (6) 

In th is work, we considered only the V passband feessell 
1990) with a zero point zpv = 3.55 x 10 -9 erg s _1 cm -2 A -1 
but the approach is absolutely general. 

To compute the SN apparent magnitude from Eq.HJone 
needs the SN spectrum at a random phase <j> to compute the 
SN flux through Eq. [5] The problem is clearly that, even in 
the best cases, SN spectral sequences have a minimum sepa- 
ration of 1 day (and even larger, particularly at later phases). 
The approach we adopted is to create a look-up table which 
enable one to construct an interpolated spectrum at any 
phase in the range spanned by the original input spectra, so 
that dangerous extrapolations are avoided. This procedure 
is described in the next section. 

2.4 Spectral Library 

The SN spectral li brary we used is co mposed by 28 spec- 
tra of SN 1998aq iBranch et al J [2 0031. merged with some 
data from SN 1994D iPatat et alJ Il99fj) and SN 1992A 
llKirshner et alJ Il993l : IStrntzeflT^ga) " spanning the phase 
range — 9 ^ t ^ 90 days from maximum light. This was pos- 
sible because these three objects are all "normal", very sim- 
ilar, Type la SNe. With the aim of increasing the generality 
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Figure 1. Dust distribution for the four spiral models considered in this work from left to right and top to bottom the total optical 
depth map is shown for models S21, S22, S31, S32 with t v (0) = 1.0 and B/T = 0.3. See the text and Table^Jfor details. 



of our approach, we devised a recipe that enables the use of 
a composite spectral library, i.e. a collection of spectra from 
different objects. One important fact to keep in mind is that 
dust extinction is a local process, i.e. it is working at the SN 
rest frame. All input spectra are first corrected to redshift 
z=0 (i.e. the SN rest frame), trimmed to the overlapping 



wavelength range and re-binned to a common wavelength 
step AA. The last two operations have the only purpose 
of making the model implementation easier. The next step 
is to ensure the consistency of the spectrophotometry with 
the LC: for each spectrum (i.e. phase) a flux normalisation 
is computed in the selected passband V, which ensures that 
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Table 1. Model parameters for the Type la SN and dust distri- 
bution used in this work. 



Component 


ro 


''max 


zo 


2-max 


la Bulge 


1.87 kpc 


8r 






la Disk 


4.00 kpc 


6ro 


0.35 kpc 


6.-0 


Dust Disk 


4.00 kpc 


6ro 


0.14 kpc 


6z 



the magnitude obtained from the spectrum (through Eqs. |S] 
and E| is equal to the one derived from the SN LC template. 

When computing the apparent flux F-p using Eq. the 
spectrum is first extinguished and then redshifted to z s i m , 
the redshift at which the simulation is focused. We note that, 
within this approach, the usable phase range is given by 
the intersection of the light curve template and the spectral 
library ones. In particular, in this work the simulation phase 
range (in the observer's rest frame) is: —5 ^ cj> ^ 80. 

2.5 The extinction coefficient C^ t (A) 

IWeineartner fc Drainel l|200lf l and lLi fc Drainel (|200lft have 
developed a carbonaceous-silicate grain model which suc- 
cessfully reproduces observed interstellar extinction, scatter- 
ing, an d infra red emission for the Milky Way. In his review, 
iDraind (l2003Tl . discussed the available evidence for dust ex- 
tinction per unit H column density in regions with different 
extinction curves. Based on t his discussion, he has renorma l- 
ized the size distributions of lWeingartner fc Drainel \2QQ ll) . 

In this work we used the extinction cross-sections per 
unit H nucleon for R Y = 3.1, which is considered to 
be appropriate for the typical diffuse Hi clouds in the Milky 
Way. The model can be used as well with the extinction 
cross-sections computed for Ry — 4.0 and Rv = 5.5. 

Further details on the grain abundance per H nucleon 
adopted for the three different Rv values and the models 
themselves can be found at Draine's web page 1 . 



3 EXPLORING THE PARAMETERS' SPACE 

The model outlined in the previous section can reproduce 
several physical situations by means of different combina- 
tions of physical and geometrical parameters. In particular: 

(i) Dust spiral arms geometry. The geometry is governed 
by the perturbation parameter uij, the number of arms N a 
and the pitch angle p. We considered five combinations to 
describe a plain exponential disk, a normal and an extreme 
spiral with two or three arms. The parameters used in each 
model and the corresponding IDs are summarized in Table 

m 

(ii) Total optical depth along the symmetry axis tv(0). 
This parameter is used to fix the total amount of dust con- 
tent of the galaxy in which SNe are being simulated. All 
the simulations in this work are for tv(0) = 1.0 because it 
is then easy to re-scale the results to different normalisa- 
tion values. M oreover this is also considered a typical value 
jBianchill2003T) . 

1 http:/ /www. princeton.edu/~draine 



Table 2. Dust geometries. 



Model 


N a 


w d 


V 


Note 


SO 









plain disk 


S21 


2 


0.4 


10° 


normal spiral 


S22 


2 


0.6 


5° 


extreme spiral 


S31 


3 


0.1 


10° 


normal spiral 


S32 


3 


0.6 


5° 


extreme spiral 




0.0 3.0 4.0 6.0 8.0 0.0 2.0 4.0 6.0 8.0 



Figure 2. Distribution of the total extinction Av for differ- 
ent galaxy inclination ang les for a S21 model (see Tab. with 
t v {0) = 1.0 and B/T = 0.3. The fraction F of objects in the in- 
clination bin considered is plotted in log scale. In each panel the 
range of inclinations considered and the total number of simulated 
SNe is shown. 



(iii) Bulge-to-Total ratio. Our assumption that the distri- 
bution of Type la SNe follows that of stars, implies that a 
fraction B/T of the SNe will be distributed in the bulge, 
while the remaining (1 — B/T) will be distributed in the 
disk. We considered three different cases for B/T: 0.0, 0.3, 
0.5. The first one describes a system in which all the simu- 
lated Type la SNe are located in the disk while the last one 
represents a system in which half of the simulated Type la 
SNe are located in the bulge and half in the disk. 

(iv) Rv- The value assumed for Rv depends on the in- 
terstellar environment intersected by the line of sight: ex- 
tinction curves which are measured in low density interstel- 
lar environments are usually characterised by small values 
of Rv (e.g. 3.1) whilst the value assumed by Rv is higher 
(4 < Rv < 6) if the line of sight is intersecting denser clouds. 
The mean extinction curve commonly used has Rv = 3.1; 
we thus adopted this value in all our simulations. 
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Table 3. Extinctions for different inclination bins. (A) indicate the mean extinction value in the V band. The subscripts / and I indicate 
that the mean refers to the whole sample and to an extinction-limited subset with A < 0.6 respectively. The table shows the figures for 
three different SN populations: B/T = 0.0 (disk), B/T = 0.3 and B/T = 1.0 (bulge). The maximum extinction value A max and the 
fraction of SNe, F, in the extinction-limited subsets are also shown. 



B/T — > 


0.0 


0.3 


1.0 


0.0 


0.3 


1.0 


0.0 


0.3 


1.0 


0.0 


0.3 


1.0 


Inclinations 


(A)f 


(A)f 


(A) f 


(A)i 


{Ah 


(A)i 


Fi 


Fi 


Fi 




^max 


^max 


0° < i < 15° 


0.15 


0.23 


0.41 


0.11 


0.13 


0.17 


0.95 


0.87 


0.68 


1.51 


1.61 


1.63 


15° i < 30° 


0.17 


0.25 


0.44 


0.12 


0.13 


0.16 


0.93 


0.85 


0.65 


1.66 


1.68 


1.71 


30° i < 45° 


0.19 


0.28 


0.49 


0.12 


0.13 


0.14 


0.91 


0.82 


0.61 


1.93 


2.04 


1.98 


45° ^ i < 60° 


0.25 


0.36 


0.61 


0.13 


0.13 


0.13 


0.86 


0.77 


0.55 


2.56 


2.56 


2.60 


60° ^ t < 75° 


0.41 


0.54 


0.87 


0.15 


0.14 


0.11 


0.77 


0.69 


0.51 


4.22 


4.28 


4.34 


75° ^ i < 90° 


1.24 


1.40 


1.78 


0.16 


0.14 


0.08 


0.56 


0.54 


0.50 


27.79 


27.34 


20.51 



3.1 Extinction distribution as a function of galaxy 
inclination 

The main geometrical parameter governing the distribution 
of the total extinction Ay is the inclination of the host 
galaxy on the sky. 

This can be clearly seen from Figure^ where the extinc- 
tion distribution is shown for six different galaxy inclination 
bins for an S21 model with -ry(O) = 1.0 and B/T = 0.3. For 
the sake of clarity we reported in Table|3the mean extinction 
for the full sample ({A)f) and for an extinction- limited sub- 
set of objects with Ay < 0.6 mag ({A)i). As expected, the 
effect of extinction in the face-on sample is modest though 
there is a significant probability, P(Av > 0.6) ~ 13% (see 
Tab.[3Jl, to have individual extinctions higher than 0.6 mag. 
In the edge-on sample the fraction of highly extinguished 
SNe is higher, as it can be seen from the extended high- 
extinction tail in the Av distribution. In this case it is 
P(Ay > 0.6) ~ 46%. These two results are easily explained 
in terms of the maximum optical depth a Type la SN can 
experience because of the combination of its position in- 
side of the galaxy and the line of sight along which the 
galaxy is seen by the observer. What is less intuitive is 
that in edge-on galaxies there is still an high probability 
that a Type la SN suffers only a small amount of extinc- 
tion P(Ay < 0.1) ~ 23%. The explanation for such effect is 
again found in geometrical arguments: the dust is confined 
on a disk which is thinner than the disk and the bulge (see 
Tab.0 on which Type la SNe are distributed. Thus, as long 
as the SN height above (or below) the galaxy plane is larger 
than the dust scale-height, the Type la SN will suffer only 
little, and in some cases even negligible, absorption. If we 
consider the extinction- limited sample (see Tab. [3] col 6), 
we see that the mean extinction still increases with the in- 
clination but only slightly: the net effect is only of a ~ 0.01 
mag increase from the face-on to the edge-on case. 

As a final test we performed the same analysis described 
above for bulge and disk SNe separately. For this purpose we 
run another two simulations, one with B/T — 0.0 to gener- 
ate all the SNe in the disk and one with B/T = 1.0 to gener- 
ate all the SNe in the bulge. We choose this approach instead 
of analysing the two SN population of the B/T — 0.3 sim- 
ulation in order to keep the analysis at the same statistical 
significance level. The figures are listed in Tab. [3 The mean 
absorption for the full sample, (A) /, increases with increas- 
ing inclination angles, as for the B/T = 0.3 case discussed 
above, but for disk SNe the mean value at each inclination is 



always smaller than for bulge SNe. The explanation for this 
apparently strange behaviour should again be found in geo- 
metrical arguments: the fact that the mean extinction value 
at any inclination is always larger for bulge SNe is due to 
the fact that those objects are located in the region where 
the dust density is higher. Indeed some of the disk SNe will 
be located at a large radial distance, where the dust density 
is very low, and thus they suffer large amounts of extinction 
only for a very small fraction of lines of sight. If we consider 
the extinction-limited subset we find that the mean extinc- 
tion value increases mildly with the inclination showing a 
0.06 mag enhancement from the face-on to the edge-on case. 
Conversely, for bulge SNe the trend is opposite: the mean 
extinction decreases with increasing inclinations. Once again 
the explanation is geometrical: the dust is confined in a disk 
which is much thinner then the bulge and thus with increas- 
ing inclination of the host galaxy the fraction of SN seen 
through the denser part of the disk is progressively reduced. 
Finally, we note that HBD98 model's predictions agree only 
partially with ours. In particular, they found that the mean 
extinction for the extinction limited sample decreases with 
increasing inclinations both for bulge and disk SNe. Given 
that their model predicts a decrease of only 0.04 mag from 
the face-on to the edge-on case for disk SNe, the difference 
is possibly due to the assumed central cut-off in their dust 
density distribution function (see also Section [2.211 . 

3.2 Extinction distribution as a function of the 
B/T ratio 

In Figure |3] we show the Av distributions for six different 
inclination bins and the two extreme values of B/T: 0.0 (thin 
line) and 0.5 (thick line) for an S21 model with tv(0) = 1.0 
and B/T = 0.3. Let us first consider, for each inclination bin, 
the peak of the corresponding distribution. With increasing 
inclination angles the relative fraction of SNe with Av < 
0.1 mag (the bin's width) in a bulge- less system and in a 
B/T = 0.5 system (F0.0/F0.5) is reverted passing from 1.24 
in face-on galaxies (0° ^ i < 15°) to 0.87 in edge-on galaxies 
(75° ^ i < 90°). This means that the presence of a bulge has 
a different effect in edge-on and face-on spirals. In the former 
case, the presence of the bulge increases the fraction of SNe 
with low extinction because the probability of a SN being 
at a larger distance from the galaxy plane is increased. On 
the contrary, in the face-on sample the presence of the bulge 
decreases the fraction of SNe with low extinction because 
the probability of a SN being above, or below, the disk (i.e. 
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Figure 3. Distribution of the total extinction Av for different 
galaxy inclination angles and the two extreme values of B/T: 
(thin line) and 0.5 (thick line) for an S21 model with Ty(0) = 
1.0 and B /T = 0.3. The fraction F of objects in the inclination 
bin considered is plotted in log scale. In each panel the range of 
inclinations considered and the total number of simulated galaxies 
is shown. For the sake of clarity, the plot has been intentionally 
limited to the range ^ Ay < 7 and to 10 -5 < F < 1. The 
number of SNe in each inclination bin is shown for both cases 
N!{B/T = 0.0), N 2 (B/T = 0.5). 



crossing a large section of the disk, where the dust is located) 
is increased. 

Another interesting feature that appears from Fig. [3 is 
that, regardless of the inclination bin considered, the B/T — 
0.5 model (thick line) appears to follow the B/T = 0.0 one 
for a while and then becomes larger at a certain value of 
Av which depends on the inclination bin considered. This 
is clearly due to the fact that, in the B/T = 0.5 model, a 
fraction of the SNe which explode in the bulge happen to 
be located also inside the disk, thus explaining the increased 
fraction of objects with higher extinction with respect to the 
bulge-less model. 



3.3 Projected radial distributions 

An interesting quantity to analyse is the projected radial 
behaviour of the mean and median absorption. The interest 
clearly arises from the possibility of comparing the predicted 
distributions with the observed ones. In Figure 3] we show 
the mean (thin line) and median (thick line) Av projected 
radial distributions for a S21 model with tv(0) = 1 and 
B/T = 0.3. The two distributions have been constructed by 
selecting all SNe located inside concentric circular bins, at 
increasing projected radial distance, and then computing, 
for each subset, the mean and the median extinction val- 
ues. As expected, the mean and median absorption values 
in the nuclear region increase with larger inclination angles 
of the host galaxy. If we consider the distribution of median 



Figure 4. Mean (thin line) and median (thick line) Av as a 
function of the projected radial distance R pr j from the centre of 
the host galaxy for an S21 model with Ty(0) = 1.0 and Rv = 3.1. 



Table 4. The effect of dust spiral geometry on SN magnitudes. 
The first column indicates the two models for which the difference 
|Atti\a| was considered. The next columns report the fraction of 
objects for which the magnitude difference obtained in the two 
models is lower than 0.025 mag, or greater than 0.2 mag and 1.0 
mag. 



Geometries 


|A] < 0.025 


|A| > 0.2 


|A| ^ 1.0 


(S0-S21) 


0.50 


0.14 


0.003 


(S0-S22) 


0.44 


0.22 


0.009 


(S0-S31) 


0.52 


0.11 


0.002 


(S0-S32) 


0.45 


0.19 


0.006 


(S21-S31) 


0.51 


0.17 


0.007 


(S22-S32) 


0.45 


0.23 


0.007 



absorption, we see that the radial functional dependency is 
quite similar for each inclination, the only feature changing 
being the peak value of the distribution. We note also that, 
in each inclination bin considered, SNe suffering high extinc- 
tion are always those projected in the very central regions. 
Finally, the fact that the median distribution is smaller than 
the mean one suggests that, at a given radial distance, there 
is a "tail" of SNe suffering high extinctions (as implied by 
the definition of mean and median). 



3.4 The spiral perturbation of the dust disk 

We performed a simple test to quantify the impact the choice 
of one dust geometry with respect to another has on the 
observed SN magnitudes. We run one Monte Carlo simu- 
lation of 10 6 SN positions and line of sights for a galaxy 
with B/T = 0.3 (the only model free parameter affecting 
the distribution of Type la SNe). We then computed the 
observed SN magnitudes for each of the five dust geometries 
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Figure 5. Effect of different dust geometries on simulated SN 
magnitudes. The three panels show the distribution of the differ- 
ence between the observed magnitude, for the same SN positions 
and line of sights, in the SO, S2 and S3 models. The fraction F 
of objects is plotted in a log. scale. The thin lines show the dis- 
tributions of AmyfSO - 521) (bottom panel), Am^fSO - S31) 
(middle panel) and Am^(S21 — S31) (top panel) whilst the thick 
lines show the distributions of Amy (SO — S22) (bottom panel), 
Am v (S0-S32) (middle panel) and Amy (S22-S32) (top panel). 
All distributions have a bin size of 0.05 mag. 



summarised in Table adopting Ty(fJ) = 1.0. In this way, 
SN positions are the same in each simulation and thus their 
magnitudes obtained in each dust geometry can be directly 
compared. In particular, we show in Figure^Jthe magnitude 
difference Amy between: the plain disk dust distribution 
SO and each of the two spiral arms ones (S21, S22; bottom 
panel); the plain disk dust distribution SO and each of the 
three spiral arms one (S31, S32; middle panel); each of the 
two spiral arms dust distributions and the corresponding 
three spiral arms ones (upper panel). From this test we can 
conclude that: 

• in each of the four cases shown in the bottom and mid- 
dle panels of Fig. the effect introduced by the spiral per- 
turbation of the dust disk with respect to the plain disk 
case on the observed magnitudes is negligible for a fraction 
of SNe which goes from a minimum of ~ 44% (when the S22 
geometry is considered), to a maximum of ~ 52% (when the 
S31 geometry is considered), being |Amv| < 0.025 mag. 

• the fraction of SNe with |Amy| ^ 0.2 mag, i.e. those 
showing a significant magnitude difference, is larger for the 
extreme spiral cases (S22, S32) than for the "normal" spiral 
cases (S21, S31). Moreover the difference between the plain 
dust disk and the spiral perturbed one is always larger in 
the two-arms cases with respect to the corresponding three- 
arms ones. The figures are shown in Table 2] Finally, the 
fraction of SNe with [Amv| > 1-0 mag is lower then 1% in 
all the cases. 



4 TYPE IA SN RATE CORRECTION FACTOR 

To exemplify how our model can be used to derive the ex- 
tinction correction factor, we will analyse a single and simple 
test case. In particular, we will consider an ideal SN search 
whose galaxy sample is made of objects of the same type, 
at the same distance and with the same properties (Rv, 
tv(0), B/T and dust geometry). This is almost equivalent 
to consider a real SN search targeted on a galaxy cluster for 
which the calculation of the correction factor is restricted to 
a single galaxy morphological class. It seems indeed plau- 
sible that for next future SN searches the statistics will no 
longer be a limitation and hence the computation of the 
SN rates in different galaxy types will become feasible. We 
thus run a simulation of 10 7 Type la SNe in a galaxy with 
tv(0) = 1.0, B/T = 0.3 at redshift z = 0.10 and with the 
dust geometry we labelled S21 (see Tab.0. 

It is obvious that the effect of host extinction is to re- 
duce the observed SN number counts and thus the derived 
SN rates will be only a lower limit if no correction is ap- 
plied. To compute the SN rate extinction correction factor 
1Z from our simulation data we just need the intrinsic mag- 
nitude distribution iVo(m) and the extinguished magnitude 
distribution N e (m). Then, 1Z will be simply given by the 
following relation: 



K : 



No(m)e(m)dm 



N e (m)e(m)dm 



(7) 



where e(m) is the SN search detection efficiency function 
(see below) and the integrals are extended over the mag- 
nitude range defined by the bright end of the magnitude 
distribution mi and the faintest magnitude for which the 
efficiency function is not null m e . Extinction corrected SN 
rates are eventually obtained simply multiplying the uncor- 
rected rates by the correction factor 1Z. The SN detection 
efficiency function, e(m), provides the sensitivity of the SN 
search as a function of the SN candidate's magnitude. The 
detection efficiency as a function of SN magnitude is usually 
a complicated function of many parameters, like the qual- 
ity of the subtracted images (seeing, transmission) and the 
technique (convolution, selection criteria) used to extract 
the SN candidates. The efficiency function is usually deter- 
mined through Monte Carlo simulations in which a number 
of artificial SNe are positioned on the se arch images (see 
iRiello et al71l2003l: ICappellaro et al.ll2004l) . These synthetic 
images are then processed with the same software used to 
perform the search on the real images and the fraction of 
recovered-over-injected artificial SNe, e = N r /N to t, is de- 
te rmined. In this work we used one of the efficiency curves 
of ICappellaro et alJ (120041) for a seeing of 0'.'90 re-scaled to 
an exposure time of 600 sec, which is a reasonable value for 
a z — 0.10 SN se arch performed with a 2m-class telescope 
(as in the case o f ICappellaro et alJ l2004). We note that the 
effect of exposure time on the efficiency curve is equivalent 
to a rigid shift of the magnitude scale. 

As an example, we have applied the method described 
at the beginning of the section. First we show in Fig. [(J 
the intrinsic, fo(m)dm (black line), and the extinguished, 
f e (m)dm (grey line) magnitude distributions normalised 
to the total number of simulated objects. The shape of 
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Figure 6. Intrinsic, fo(m)dm (black line), and extinguished, 
f e (m)dm (grey line), magnitude distributions. The connected 
solid dots represent the SN detection efficiency curve used in this 
work (see the text for further details). 

fo(m)dm is linked to that of the SN LC: at phases cj> > 
30 — 40d the LC is much less steep than at earlier phases 
and thus, for a given magnitude bin, there is a wide range 
in phase (i.e a large number of SNe) contributing to that 
bin. The spike and wiggles at 21 mv ^ 22 are mainly due 
to the fact that in this simulation we used only one specific 
Type la LC. We expect them to smooth out when Type la 
SN diversity is taken into account by combining simulation 
with different SN properties. From our calculation we find 
72 = 1.27, i.e. ~ 27% of SNe are not detected by the search 
because of host galaxy extinction. Given the impressive size 
of the tail of the f e (m)dm in Fig. |S| it is interesting to 
look at the intrinsic fraction (without taking into account 
the efficiency curve) of SNe that are extinguished beyond 
the fainter magnitude reached in the absence of extinction 
(mv — 22.08 in the case shown in Fig. [SJ: the tail includes 
~ 24% of the simulated SNe. 



5 DISCUSSION 

The recipe described above to derive the correction factor 
72 clearly depends on the model parameters. To assess their 
importance on the final result we performed several simula- 
tions to probe its dependence from: (1) dust geometry, (2) 
size of the bulge, (3) dust content and (4) dust type. 

1. We computed the correction factor for each dust ge- 
ometry described in Tab. fusing Ty(0) = 1.0, Rv = 3.1 and 
B/T = 0.3. We found that 72 = 1.27 in all cases. This result 
demonstrates that the distribution of dust does not play a 
major role and thus it is possible to adopt the simplest dis- 
tribution, i.e. the plain double exponential disk SO. Though 
this is true when the extinction correction TZ is concerned, 
we note that a detailed dust distribution should be included 



when studying the distribution of extinction and its effect on 
SN properties/distributions (see Sec.[3Jl. This result is not 
unexpected because (a) the spiral perturbation of the dust 
distribution creates regions of enhanced (the spiral arms) 
and lower density with respect to the plain exponential disk; 
and (b) Type la SNe are distributed (in the simulations) on 
a plain exponential disk, i.e. without a spiral perturbation. 
The effect of dust over- and under-densities therefore can- 
cels out. Conversely we expect that the dust distribution 
will play a more significant role for CC SNe because they 
are also more concentrated in the spiral arms. 

2. We computed the correction factor for 0.0 ^ B/T ^ 
0.5 using dust model S21 with ry(O) = 1.0 and Rv = 3.1. 
We found that the corresponding range spanned by the ex- 
tinction correction factor is 1.22 ^ TZ ^ 1.31. Hence, the 
fraction of Type la SNe that are located in the bulge, which 
is governed by the B/T parameter, plays a role not only 
on the extinction distribution (see Sec. 13.21 but also on the 
extinction correction factor 72. 

3. We computed the correction factor TZ(tv) for dif- 
ferent values of the total face-on optical depth Ty(0) using 
dust model SO with Ry = 3.1 and B/T = 0.3. We found: 
11(0.5) = 1.16, TZ(1) = 1.27, 72,(2.5) = 1.57, 72(5) = 1.91 
and 72(10) = 2.35. 

4. We computed the correction factor for Rv ~ 3.1, 
4.0 and 5.5 using dust model SO with tv(0) = 1.0 and 
B/T = 0.3. We recall that these are the values for which the 
extinction coefficients C^ t (A) are available (see Sec. 12.51 . 
We found that the corresponding values of the extinction 
correction factor are 72 = 1.27, 1.31, 1.34. 

When the determination of 72 is concerned, although 
the detailed geometry of dust can be neglected, the dust 
content, properties and the bulge size still play a major role. 
The computation of 72 for a real SN search requires indeed 
several simulations: one for each "galaxy class" which is de- 
fined by the tuple (B/T, Ry,Tv). If the SN rates cannot be 
computed (and hence corrected) for each galaxy class sepa- 
rately, a global correction factor (72) should be derived from 
a weighted mean of the 72i obtained for each galaxy class i, 
with weights given by the fraction of objects in each galaxy 
class. 



6 SUMMARY 

We presented a new Monte Carlo approach to the problem of 
determining the effect of host galaxy extinction on observed 
SN properties and to compute the extinction correction fac- 
tor for SN rate. Although we provided a detailed discussion 
of the role played by the model free-parameters, our main 
purpose was to outline the recipe to derive the extinction 
correction factor 72 for Type la SN. Our main results are 
that (1) the value of 72 is almost insensitive to the partic- 
ular spiral geometry of dust and hence a plain exponential 
disk is a good approximation. (2) the value of 72 strongly 
depends on the galaxy physical properties, in particular the 
total face-on optical depth ry(0), the value of Rv and the 
size of the bulge B/T. If we adopt a Milky Way-like dust 
(Rv = 3.1), B/T = 0.3 and the range of most likely total 
face-on optical depths 1 tv(0) «S 5 jBianchill200a) . the 
corresponding range spanned by the extinction correction 
factor is 1.27 < 72 «S 1.91. 
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The computation of 1Z for a real SN search involves 
a number of complications due to the heterogeneity of the 
galaxy sample, the redshift dependence of galaxy properties 
and finally the variety of Type la SNe. In particular, this 
work did not take into account neither the intrinsic spread 
in Type la SNe absolute magnitudes nor the correlation be- 
tween the light curve shape (characterised, e.g. by the Amis 
parameter) and the morphological types of SN hosts. Such 
simplifications were possible only because the aim of our 
work was to provide an overview of the method rather than 
discussing the application to a real case, which will be in- 
stead the subject of a future work. This application clearly 
requires a detailed study of the observed galaxy sample and 
the definition of the corresponding model free-parameters 
which was beyond the objective of the present study. Fi- 
nally, we note that, though the approach described in this 
work can be used also for core collapse SNe, its feasibility is 
limited both by the lack of adequate spectral libraries and 
by the extreme heterogeneity of these objects. 
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